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Chapter  VIII 


Three-Dimensional  Flames 


The  Flame  as  a  Hydrodynamic  Discontinuity. 

The  plane  premixed  flame  discussed  in  Chapters  II  &  III  is  an  idealization 
seldom  approximated,  since  in  practice  the  flame  is  usually  curved.  A  bunsen 
flame  is  quite  unavoidably  so,  but  even  under  circumstances  carefully  chosen 
to  nurture  a  plane  state,  instabilities  can  lead  to  a  three-dimensional  structur 
Such  flames  have  been  extensively  studied  (Markstein  1964,  p.  7)  using  what  may 
be  called  a  hydrodynamical  approach,  a  brief  description  of  which  will  provide 
an  appropriate  introduction  to  our  subject. 

On  a  scale  that  is  large  compared  to  its  nominal  thickness  X/c  M,  the 

P 

flame  is  simply  a  sheet  across  which  there  is  a  Jump  in  temperature  and  density 
subject  to  Charles'  Law,  as  is  appropriate  for  an  isobaric  process.  Deformation 
of  the  sheet  from  a  plane  is  associated  with  pressure  variations,  of  the  order 
of  the  square  of  the  Mach  numbers  (cf.  Sec.  1.5),  in  the  hydrodynamic  fields. 
These  small  pressures  jump  across  the  sheet  in  order  to  conserve  normal  momentum 
flux.  Because  Euler's  equation  for  small  Mach  number  hold  outside  the  flame 
(cf.  end  of  Sec.  3),  the  temperature  and  density  do  not  change  along  particle 
paths;  so  that  for  a  flame  travelling  into  a  uniform  gas  the  temperature  and 
density  ahead  are  constant  and  the  flow  is  irrotational.  The  flow  behind  is 
stratified,  however,  since  flame  curvature  generates  both  vorticiy  and  non- 
uniform  temperature  jumps.  Nevertheless  variations  in  temperature  from  the 
adiabatic  flame  temperature  are  usually  neglected  everywhere,  a  matter  to 
which  we  shall  return  later.  On  the  other  hand,  vorticity  generation  cannot 
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be  ignored  so  that,  although  Euler's  equations  hold  behind  the  flame,  the 
flow  is  not  potential  there. 

The  flow  fields  on  the  two  sides  of  the  flame  are  coupled  by  conservation 
of  mass  and  momentum  fluxes  through  the  sheet.  Indeed,  if  V  is  the  speed  of 
the  discontinuity  back  along  its  normal  n  (Fig.  l)  we  have 


d,2,3)  p1(vnl 


where  vq  =  yn  and  y^  =  y  -  vnp.  The  flame  speed  is  defined  to  be  the 
normal  gas  speed  immediately  ahead,  as  measured  in  a  local  frame  fixed  at  the 
discontinuity,  i.e.  +  V.  If  this  is  specified  the  hydrodynamic  problem 

can  in  principle  be  solved  and  the  locus  of  the  flame  determined.  [Equations 
(l),  (2),  (3)  become  four  scalar  equations  for  v2>p2  once  p2  i‘s  known» 
and  that  follows  from  the  flame  temperature  Tg. 

Evaluation  of  the  flame  speed  from  a  combustion  analysis  of  the  sheet 
has  often  been  sidestepped,  and  replaced  by  one  of  several  hypotheses.  The 
simplest  is  that  the  speed  is  constant.  In  this  way  Landau  (19M»)  considered 
infinitesimal  disturbances  of  a  plane  flame  and  concluded  it  is  always  unstable, 
a  result  that  will  be  examined  in  Ch.  XI. 

The  already  nonlinear  flow  field  is  only  complicated  by  the  combustion, 
even  for  a  constant-speed  flame.  Accordingly  an  additional  simplication, 
often  adopted,  is  that  the  flow  is  not  affected  by  the  flame.  This  can  be 
justified  only  when  the  heat  released  by  the  combustion  is  small  compared 
to  the  thermal  energy  of  the  fresh  mixture  (Sec.  1.5),  so  that  the  temperature 
Jump  across  the  sheet' is  negligible.  Such  is  not  a  common  characteristic  of 
actual  flames,  but  the  approach  can  provide  insight  into  their  qualitative 


nature . 
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Consider,  for  example,  a  stationary  flame  located  in  the  shear  flow  of 
Fig.  2,  representative  of  that  from  a  bunsen  tube.  Since  vn  is  equal  to  the 
flame  speed,  assumed  constant,  the  arigle  a  decreases  as  u  increases  and 
there  is  a  cusp  on  the  centerline.  The  resulting  flame  shape  has  some  resemblance 
to  that  of  the  inner  cone  of  a  bunsen  flame.  (Even  for  a  parabolic  velocity 
distribution,  a  good  approximation  to  the  flow  from  the  tube  when  that  is  unaffected 
by  the  combustion,  the  shape  of  the  flame  can  only  be  obtained  numerically. ) 

In  practice  the  tip  is  usually  rounded  so  that  the  flame  speed  at  the 

centerline  is  equal  to  v  and  must  have  a  smaller  value  elsewhere.  (A 

max 

detailed  study  of  flame  tips  is  presented  in  Ch.  IX. )  Observations  of  this 
nature  and  dissatisfaction  with  the  stability  prediction  of  Landau  have 
encouraged  more  sophisticated  hypotheses,  in  particular  that  of  Markstein 
discussed  at  the  end  of  Sec.  III. 6.  We  shall  not  pursue  these  here  since 
rational  analysis  of  the  burning  zone  (which  is  embedded  within  the  hydro- 
dynamic  discintinuity)  determines  the  actual  flame  speed  from  first  principles. 


2.  Slow  Variation  and  Near  Equidif fusion. 


As  discussed  in  Sec.  II. 7,  the  combustion  of  a  two-component  mixture  that 


is  well  removed  from  stoichiometry  can  be  characterized  by  a  single  mass 
fraction  Y  related  to  the  component  which  is  consumed  at  the  sheet.  To  be 


sure,  the  argument  was  given  for  plane  flames,  equal  molecular  masses,  and 
unit  Lewis  numbers  and  stoichiometric  coefficients;  but  all  these  restrictions 
are  readily  seen  to  be  inessential  (as  also  is  the  number  of  components).  The 
single  Lewis  number  is  then  properly  interpreted  as  that  of  the  deficient 


component  -  oxidant  in  a  rich  flame,  fuel  for  a  lean  one. 
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Appropriate  equations  for  our  discussion  of  three-dimensional  unsteady 
motion  are  therefore 


(U,5) 

(6) 

(7) 

(8) 


pT  =  1,  3p/3t  +V  *(pv)  =  0, 


p(8Y/3t  +  v-V  Y)  -2._172Y  =  -fi. 


p 

p(  3v/  3t  +  v'Vv)  =  -  7p  +  lP[  7  v  +v(V*v)/3], 


p(  3T/3t  +  v-  VT)  -V  T  =  n 


where 

(9) 


n  •  AYe 


-e/T 


These  have  been  rendered  dimensionless  as  in  Sec.  1.6  and  the  undisturbed 
pressure  has  been  taken  constant  (i.e.  pc  =  l),  with  p  the  variations  from  it. 

On  the  length  scale  X/c^M  now  being  used  the  flame  is  no  longer  a 
sheet  but,  for  large  activation  energy,  the  reaction  zone  within  it  can  be 
considered  so.  We  now  take  the  x-axis  instantaneously  parallel  to  the  normal 
to  this  sheet  (and  therefore  to  the  hydrodynamic  discontinuity)  at  the  point 
of  interest  and  introduce  the  new  variable 


(10) 


n  =  x  -  x#(t) , 


as  was  done  in  the  discussion  of  plane  deflagration  waves  in  Sec.  III. 3.  As 
usual  x#  denotes  the  position  of  the  sheet.  Everywhere  in  the  above 
equations  3/3t  and  3/3x  should  now  be  replaced  by 


(11) 


3/3t  -  x„3/3n  and  3/3n, 


respectively. 
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to  0+  to  obtain 


<12>  I^PliS-  *  <\  -  *  V^!E  d”  -  §  I  Of  +  *  T> 

where  as  before  the  subscript  j_  denotes  the  component  perpendicular  to  n, 
i.e.  in  the  flame  sheet,  and  we  have  used  the  fact  that  Y  vanishes  behind 
the  sheet.  The  middle  term  on  the  left-hand  side  can  be  integrated  by  parts 


0+„2,,-l 


In  analysing  the  resulting  equations  for  the  limit  Q-*»  some  thought 
must  be  given  to  the  flame  temperature,  i.e.  the  temperature  immediately 
behind  the  sheet.  It  is  clear  from  the  results  of  Chs.  II  &  III  (cf.  Sec.  III.l) 
that  there  must  be  a  general  correlation  between  this  temperature  and  the 
flame  speed  in  the  sense  that  0( 1/0 )  changes  in  the  former  are  associated 
with  0(l)  changes  in  the  latter.  We  may  therefore  expect  that  if  the  flame 
temperature  varies  by  larger  amounts,  either  temporally  or  spatially,  then 
the  flame  will  vary  in  an  extreme  fashion  through  large  acceleration  or 
curvature;  and  such  behavior  will  introduce  serious  difficulties  into  the 
analysis.  Only  in  the  discussion  of  explosions  in  Ch.  XII  will  these 
difficulties  be  confronted,  and  then  only  in  relatively  simple  circumstances. 
Here,  as  in  Ch.  Ill,  we  shall  restrict  ourselves  to  physical  situations  where 
the  temperature  of  the  burnt  mixture  in  the  neighborhood  of  the  flame  is 
everywhere,  and  at  all  times,  within  0(l/e)  of  the  adiabatic  flame 
temperature.  Necessary  restrictions  are  established  by  the  following  argument 
(Buckmaster  1979a). 

As  in  Sec.  III. 3,  we  can  calculate  the  flame  temperature  from  the  overall 
change  in  partial  enthalpy  H  =  Y  +  T  in  the  flame  up  to  and  including 
the  reaction  zone.  Thus  equations  (6)  and  (8),  with  the  modifications  (ll), 
are  added  to  eliminate  the  reaction  terms  and  then  integrated  from  n  =  -» 


to  yield 
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pH(vn  "  V\JTL  h  to  Ep(vn  ‘  *#)]dn> 


in  which 


3n 


[p(vn  -  *,)]  -  -  |£  -v^fpy^ ,  p(vn  -  i.)|M  *  P1<Tni-i.)-,^r||  «jj(pyj_! 


according  to  the  equation  of  continuity  (5).  Equation  (12)  may  therefore  he 
rewritten 


(13)  Pl(T  -i,)(T.  -  H,)  -  f  I  .  *  T)  +  I.<f*  + 

X  Or 


+  ? i  ’  PY | )  -  -jr(pH)  *  (pHv. )]  dn 


‘1  Ki_L  at 


IL  ~J_ 


where,  as  usual,  T#  denotes  the  temperature  at  the  flame  sheet. 

By  assumption  the  left-hand  side  of  equation  (13)  is  at  most  0(l/e) 


so  that  the  right  side  must  he  similarly  hounded.  The  only  way  to  guarantee 


this  within  a  general  framework  is  to  make  the  terms  in  and  3/3t 


separately  small,  which  has  heen  accomplished  in  two  ways.  We  confine 
attention  to  disturbances  which  vary  over  times  and  distances  of  order  0  . 
Such  slowly  varying  flames  are  a  generalization  of  those  introduced  in 
Sec.  III. 5*  Alternatively,  we  restrict  the  Lewis  number  to  values  close  to 
one ,  i . e .  write 


dM 


X  1  *  1  -  X/e  with  X  =  0(1), 


insist  that  the  state  1  does  not  vary  along  the  flame  and  impose  the  integral 
(15)  H  =  H1  +  0(1/0) 


which  is  then  a  possibility.  Such  near-equidiffusional  flames  are  a  generalizati 
of  those  introduced  in  Sec.  III. 3. 
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The  full  generality  of  the  above  analysis  has  not  appeared  elsewhere 
although  its  essence,  for  shear  flows,  is  contained  in  Buckmaster  (1979-a). 

The  rest  of  the  chapter  is  concerned  with  the  detailed  formulation  of  each 
of  these  two  kinds  of  flame. 

3.  The  Basic  Equation  of  Slowly  Varying  Flames. 

Slowly  varying  flames  can  be  characterized  by  the  scaling 
(16)  (y,z,t)  =  e(n, s,t), 

but  it  should  be  recognized  from  the  outset  that  the  actual  lengths  involved 
are  quite  moderate.  In  Sec.  II. 4  we  saw  that  the  preheat  zone  typically  has 
a  thickness  of  less  than  0.1  mm.  So  that,  even  for  6  as  large  as  100, 
appreciable  changes  occur  in  less  than  1  cm. 

Flame-sheet  disturbances  on  the  8-scale  will  be  associated  with  field 
disturbances  at  distances  n  =  0(e),  where  the ‘governing  equations  are 
fundamentally  different  from  those  where  n  =  0(l).  The  limit  0  » 

therefore  generates  a  structure  with  at  least  three  zones:  the  flame  sheet 
described  on  the  scale  n  =  0(l/e),  to  which  all  the  reaction  is  confined; 
a  diffusion  zone  described  on  a  scale  n  =  0(l),  called  a  preheat  zone  in 
front  of  the  flame;  and  an  ideal-fluid  region  described  on  the  scale  n  =  0(9). 
Description  of  the  latter  is  the  hydrodynamic  problem  discussed  in  Sec.  1. 

The  stratification  behind  the  flame  which  is  an  essential  feature  of  this 
description  must  eventually  be  smoothed  out  by  diffusion,  so  that  additional 
structure  is  described  in  a  scale  much  larger  than  n  =  0(0);  but  that  will 


not  be  considered. 
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Anticipation  of  the  matching  conditions  imposed  by  the  ideal-fluid 
region  enables  the  diffusion  zone  to  be  analyzed.  These  conditions  are 


(17)  n  -*•  +»:  (p,T,Y,y,p)  ->  (p1,T1,Y1>y1,p1); 

(18)  n  -*■  exponential  growth  not  permitted. 


The  state  1  is  that  immediately  ahead  of  the  hydrodynamic  discontinuity,  and 

is  to  be  considered  a  known  function  of  n»S,T.  Our  aim  is  the  same  as  that 

in  Sec.  III. 5,  namely  to  calculate  various  quantities  in  the  overall  enthalpy 

balance  (13)  and  hence  derive  an  expression,  in  terms  of  the  state  1,  for  the 

perturbation  of  the  flame  temperature.  (The  change  of  units  introduced  in 

Sec.  III. 3  is  not  appropriate  here  since  p, ,  which  replaces  p  there, 

—00 

is  in  general  variable;  hence  the  differences  in  some  of  the  subsequent  formulas.) 
To  leading  order  the  governing  equations  become 


(19,  20) 

(21) 


pT  =  1,  9[p(V  +  v  )]/3n  =  0, 

n 


p(V  +  vn)3Y/3n  -£-132Y/3n2  =  -Q  , 


(22,  23)  p(V  +  vn)3vn/3n  =  -3p/3n  +  (Uf/3)3  vQ/3n  ,p(V+vn)3v^/3n  =^3  v^/Sn  , 


(21*) 


p(V  +  v  )  3T/3n  -  3  T/3n^  =  a 
n 


where  -x*  has  now  been  identified  with  V,  the  speed  of  the  hydrodynamic 
discontinuity.  Apart  from  notation,  equations  (19),  (20),  (2l)  and  (24) 
are  identical  to  those  for  a  plane  flame  so  that  we  may  immediately  write 
(cf.  Sec.  III. 5) 


(25) 


p(V  +  vn)  =  p1(V  +  vn])  =  Mn  (say), 


] 
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£M  n  M  n  M  n  -1 

(26)  Y  =  Y]_(l  -  e  n  ),  T  =  ^  +  Yxe  n  ,  p  =  +  Y^e  n  )  for  n  <  0, 

(27)  Y  =  0,  T  =  +  9  ^’t£o(ri1  ?,t)  ,p  =  for  n  >  0. 

The  mass  flux  may  also  he  called  the  burning  rate.  Expressions  for  v 

and  p  can  also  be  written,  but  we  shall  not  need  them.  On  the  other  hand,  v 
is  now  needed;  since  it  does  not  .Jump  across  the  hydrodynamic  discontinuity,  it 
can  be  expected  to  stay  constant  over  the  whole  range  of  n: 


(28) 


Formal  proof  comes  from  equation  (23),  which  holds  even  in  the  flame  sheet 
and  possesses  only  constant  vectors  as  bounded  solutions. 

The  perturbation  t^  in  the  flame  temperature  may  now  be  calculated 
from  the  enthalpy  balance  (13)  since  all  other  quantities  have  been 
determined,  in  terms  of  V,  to  sufficient  accuracy.  Hie  result  is  quite 
complicated,  involving  terms  in  d/dt  and  of  p^T^.Y^  and  v^. 
Considerable  simpliciation  occurs  when  the  upstream  temperature  and  composition 
of  the  mixture  (on  the  hydrodynamic  scale)  are  uniform,  since  then  p^, 
and  Y^  are  constant,  as  we  shall  see  at  the  end  of  this  section.  Only 
the  simpler  result  will  be  written  here. 

The  0(9  1)  terms  in  the  integrand  ( 13 )  are  now 

-  H)]  +  ^[pv^T*  -  H)] 


=  {[hi 


111 


■?_L 


M  ]— 
nJM 
n 


3n 


V  •  v  1 

~1  Ul 


-f  M  n  M  n 

v  (  n  n  \ 

Yj(e  -  e  ) 


Tl  + 


Yle 


M  n 
n 


which  may  be  replaced  by 
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{- 


~[M 
M  n 
n 


*  MJ  * 


!j_  ’V 


£M  n 
n 


Yl(e 


M  n 

n  \ 
-  e  ) 


T1  +  Yle 


M  n 
n 


under  integration  by  parts.  Here  the  dot  denotes  differentiation  with  respect 
to  x  and  is  the  gradient  operator  with  respect  to  n, C  .  For  Vj^  =  0 

the  calculation  is  the  same  as  for  a  plane  flame,  so  that  it  is  not  surprising 
to  find  the  generalization 


(29) 


t  =  -bM 


,-3 


[M 


5rJU'-±M.> 


-  MnVr  Yjjl  where 


b  =  Y. 


J1  a 


l  o 


VY1° 


do 


of  the  result  (ill. 39)  from  the  integration  (13).  (A  factor  is  now 

missing  from  b  because  of  the  different  unit  for  in  and  the  use 

of  1/6  as  expansion  parameter. ) 

Another  relation  between  the  burning  rate  M  and  the  perturbation  t 

n  oo 

of  the  flame  temperature  follows  from  the  structure  of  the  flame  sheet, 
which  is  described  on  the  scale  n  =  0(1/6).  The  analysis  is  the  same  as  in 
Ch.  3  since,  to  leading  order,  the  equations  for  p,  Y  and  T  are  those  for 
a  plane  flame.  We  conclude  that 
W2T* 

(30)  M  =  e  with  T*  =  H.  , 

n  *  x 


cf.  equation  (ill. 26),  if  the  mass  flux  M  on  which  units  have  been  based 
is  chosen  to  be  the  burning  rate  of  a  plane  flame  moving  into  the  state 


(31)  M  =  /2XD  T*  exp(-9/2T*)/Y19, 

cf.  equation  (lll.l6).  The  two  results  (29),  (30)  now  give  us  the  basic 


equation 
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(32)  b[fon  +  Yil-Vj_  Mn  -  Mnyj_  -yj  +  2T^  te  Mn  -  0 

of  slowly  varying  flames,  which  should  he  compared  with  its  plane  from  ( 111.1+0) 
The  equation  is  due  to  Buckmaster  (1977),  whose  derivation  has  been  simplified. 
Restricted  forms  of  it  have  been  given  by  Sivashinsky(l97** ,  1975)  and,  in  a 
linearized  context,  by  Eckhaus  (l96l).  Interpretation  involves  the  concept 
of  flame  stretch  to  be  introduced  in  the  next  section.  V7e  shall  finish  the 
present  section  with  a  brief  discussion  of  the  ideal-fluid  region  where  n  =  0( 
To  the  slow  variables  n,C,T  we  now  add 

(33)  £  =  x/6. 

Then  p(=  1/T),  T  and  Y  are  found  to  be  constant  on  particle  lines,  to 
leading  order,  while  the  velocity  and  pressure  fields  satisfy 

(34)  V  •  y  ■  0,  p( 3y/ 3t  +  v  -Vv)  =  -  Jp, 

where  V  is  now  the  gradient  operator  with  respect  to  £,n,£.  These  are 
Euler's  equations  for  an  incompressible  (but  not  necessarily  constant-density) 
ideal  fluid,  i.e.  one  devoid  of  viscosity  and  heat  conduction.  Particles 
retain  the  values  p^,T^  and  Y^  they  acquire  far  upstream  until  they 
reach  the  hydrodynamic  discontinuity,  where  they  exchange  them  for  the  values 

(35)  p2  =  H"1,  T2  =  Y2  =  0, 

which  they  then  carry  downstream.  The  ideal-fluid  region  is  coupled  to  the 
flame  through  the  velocity  and  pressure, in  a  manner  that  will  now  be  seen  in 
detail  for  uniform  temperature  and  composition  upstream. 
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p,  T  and  Y  are  all  constant  on  either  side  of  the  hydrodynamic 
discontinuity.  In  particular  there  is  no  thermal  stratification  on  the  down¬ 
stream  side,  in  accordance  with  a  common  assumption  mentioned  in  Sec.  1.  The 
problem  is  to  find  v  and  p  subject  to  the  jump  conditions  Cl),  (2),  (3) 
and  the  evolution  equation  (32),  which  may  be  thought  of  as  five  equations 
for  v^,  pg  and  V  in  terms  of  the  state  1.  Clearly  its  complete  analysis 
would  be  formidable. 

4.  Flame  Stretch. 

In  view  of  the  difficulties  in  giving  a  complete  analysis  of  a  three- 
dimensional  flame  (even  when  it  is  slowly  varying),  it  is  natural  to  try  to 
identify  special  characteristics  that  play  particularly  important  roles  in 
the  understanding  of  flame  behavior.  Flame  speed  and  temperature ‘are 
elementary  examples  of  such  characteristics  that  have  already  been  identified. 

A  more  subtle  characteristic,  attributed  to  Karlovitz  (1953),  is  flame  stretch. 

In  order  to  define  this  in  an  unambiguous  fashion  it  is  first  necessary 
to  define  a  flame  surface,  i.e.  a  sheet  that  characterizes  the  location  of  the 
reaction.  For  large  activation  energy  the  reaction  zone  is  such  a. surface 
when  viewed  on  the  scale  of  the  preheat  zone.  Lewis  &  von  Elbe  (1961,  p.  225) 
use  the  locus  of  inflection  points  of  the  temperature  field  near  the  burning 
zone,  a  definition  that  coincides  with  the  preceding  one  in  the  limit  0  -»•  w>  . 
If  the  flame  can  be  viewed  as  a  hydrodynamic  discontinuity,  as  in  Sec.  1,  then 
the  discontinuity  itself  is  a  flame  surface.  In  all  these  examples  and  any 
other  for  which  the  present  discussion  is  meaningful,  a  flow  velocity  is 
defined  on  each  side  of  the  surface  such  that  v  is  continuous  across  the 


surface. 
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Consider  a  point  on  the  flame  surface  which  moves  in  it  with  the  velocity 
.  A  set  of  such  points  forming  an  infinitesimal  surface  element  of  area 
S  will,  in  general,  deform  during  the  motion,  so  that  S  will  vary  with  time. 
If  S  increases  the  flame  is  said  to  be  stretched,  whereas  if  S  decreases 
the  flame  experiences  negative  stretch  and  is  said  to  be  compressed.  A  precise 
measure  of  the  stretch  is  the  Karlovitz  number 


(36) 


_  1  dS 
"  S  dt  * 


Certain  configurations  for  which  stretch  is  an  important  characteristic  will 
be  discussed  in  detail  in  Ch.  X.  Here  we  small  give  two  elementary  examples. 

Consider  a  stationary  flame  surface  which,  for  simplicity,  is  locally  plane, 
inclined  at  an  angle  a  to  the  horizontal  (Fig.  3).  If  the  flame  is  located  in 
a  horizontal  shear  flow  U(y),  then 

(37)  ~  “■  =  sina  cosa  U'(y). 


Note  that  if  a  =  0  or.  -rr/2  there  is  no  stretch  because  then  the  tangential 
velocity  does  not  change  along  the  surface.  The  second  example  is  an  expanding 
spherical  flame.  If  the  radius  is  R(t),  then 


(38) 


I  as  _  2  • 

S  dt  ■  R  " 


Stretch  of  this  type  is  discussed  under  somewhat  more  general  circumstances 
in  connection  with  the  quasi-plane  flames  of  Sec.  III. 8. 

In  these  two  examples  the  stretch  could  be  determined  by  inspection, 
but  we  new  need  a  general  expression  for  it.  That  comes  from  connecting  it 
with  stretch  in  a  plane.  Thus,  as  an  element  of  the  flame  surface  passes 
through  a  designated  point  it  experiences  the  same  stretch  as  its  projection 


on  the  tangent  plane  at  that  point,  which  is  moving  with  velocity  Vj_.  The 
plane  element,  however,  has  stretch 


(39) 


1  dS 
S  dt 


as  is  easily  seen  from  the  divergence  theorem,  so  that  this  formula  holds 
also  for  the  surface  element. 

In  the  first  example  above  the  flame  surface  is  its  own  tangent  plane  at 
every  point  and  we  have 


Yj^=  (U  cos  a,  U  sina  cos  a,  0), 


from  which  the  result  (37)  follows.  In  the  second  example 


Yj_=  (0,  Ry/R,  Rz/R), 


if  the  x-axis  is  taken  along  the  normal  at  the  point  considered,  and  the 
result  (38)  follows. 

Stretch  is  a  kinematical  concept.  For  a  slowly  varying  flame  with  uniform 
temperature  upstream  it  may  be  related  to  the  flame  speed.  The  burning  rate  !■' 
is  then  essentially  the  flame  speed  v  ^  +  V,  differing  from  it  only  by  a 
constant  factor  p^.  The  basic  equation  (32),  which  now  reads 

(l,0)  I  f  -  ih  »  V 

n 

therefore  expresses  the  stretch  in  terms  of  the  flame  speed.  Note  that  stretch 
is  the  same  for  the  hydrodynamic  discontinuity  as  for  the  reaction  zone,  but 
flame  speed  relates  to  the  former  and  not  the  latter. 

Flame  thickness  is  another  concept  that  can  usefully  characterize  a  flame. 
Although  a  plane  deflagration  wave  is  not  confined  to  a  zone  of  finite  thickness. 


to  say  that  it  has  a  thickness  X/c  M  obviously  has  meaning.  The  implication 

P 

is  that  while  an  infinite  distance  is  required  for  the  temperature  to  rise 
from  the  cold  value  to  the  hot,  most  of  the  change  takes  place  over  a  small 
multiple  of  X/c^M. 

For  a  general  flame,  it  is  natural  to  define  the  local  thickness  as 
X/c  where  is  the  normal  mass  flux.  Indeed,  for  the  slowly  varying 

flame  the  dimensionless  Mn  is  the  decay  coefficient  for  T  in  the  preheat 
zone,  as  the  second  of  formulas  (26)  shows,  so  that 


(1*1) 


£  =  M 


-1 


n 


is  properly  called  the  non-dimensional  thickness.  If  now  the  elemental  volume 
(42)  A  =  £S 


generated  by  S  is  introduced,  the  basic  equation  becomes 

(43)  7-  ~  =  2T?b-1M^  £n  M  , 

A  dr  *  n  n’ 

expressing  what  may  be  called  voluminal  stretch  in  terms  of  flame  speed  or 
thickness.  The  results  (4o)  and  (1*3)  were  originally  derived  by  Buckuacter  (l 
Note  that  the  sign  of  b  depends  on  the  Lewis  number  ji,  ,  being  positive 
for -/.<  1  and  negative  for  ~j>  1.  When  /  =  1  +  o(  1/© )  so  that  b  =  0(l/e), 
equation  (43)  suggests  that  changes  occur  on  the  scale  of  t  rather  than  t  , 
although  it  then  fails  to  provide  a  valid  description.  In  Sec.  8  we  shall 
develop  the  theory  of  near-equidiffusional  flames,  and  fill  the  void. 


We  shall  end  the  section  with  a  few  remarks  about  the  contemparary 
significance  of  flame  stretch.  Although  the  concept  leads  to  an  elegant 
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description  of  slowly  varying  flames  (especially  when  extended  to  vo luminal 
stretch,  about  which  the  literature  is  silent),  the  importance  attached  to  it 
derives  from  its  far-ranging  use  as  an  intuitive  tool  in  the  prediction  and 
explanation  of  flame  behavior  (Lewis  &  von  Elbe,  19&1 ).  This  is  not  the 
appropriate  place  to  discuss  the  arguments  used,  which  are  in  fact  no  more 
than  speculation.  We  Just  mention  here  the  fundamental  claim  (sometimes  modified 
in  an  inconsistent  fashion  to  fit  the  facts)  that  stretching  a  flame  causes  it 
to  decelerate.  The  mathematical  evidence  for  such  a  simple  picture  is  mixed. 
Certainly  the  present  results  for  slowly  varying  flames  are  not  favorable: 
equation  (32)  shows  that  deceleration  is  associated  with  negative  stretch.  But 
there  is  supportive  evidence  from  near-equidif fusicnal  flames  in  certain  special 
circumstances.  These  will  be  discussed  in  Ch.  X. 

5.  Ncar-Equidiffusional  Flames. 

The  nonun.ifonuity  as  ■)'  -*  1  revealed  by  equation  (32)  makes  it  clear  that 
a  different  formulation,  distinct  from  that  for  slowly  varying  flames,  is 
needed  when  heat  and  reactant  diffuse  at  comparable  rates.  The  essential 
characteristics  of  such  a  formulation  have  been  identified  in  Sec.  2,  in 
particular  equations  (l4)  and  (15);  we  shall  now  present  a  detailed  mathematical 
treatment.  The  starting  point  is  the  system  of  equations  (3-9)  with 

(4*0  A  ~  (Y^G2/2tJ|)  exp(0/T*)  with  T*  =  H  , 

corresponding  to  the  choice  (3l),  with  '/ =  1,  for  M.  Recall  that  the  state 
1  does  not  vary  along  the  flame. 


t 


We  seek  en  asymptotic  development,  in  which  spatial  and  temporal  variations 
in  the  flame-sheet  location  are  0(l).  As  discussed  in  Sec.  2,  limitation  of 


-17- 


the  allowable  variation  in  flame  temperature  then  imposes  the  restriction 
(l4)  and  leads  to  the  integral  (15),  which  we  now  write 

(45)  H  =  T*  +  e_1h  +  o(e"2) 

outside  the  reaction  zone.  Equations  (6),  (8)  may  then  be  written 
(46,47)  p(~  +  y?)T  =  V2T,  p(-^  +  Y  -  J  )h  =  V2h  +  XV2T 

there,  where  p,  T,  v  now  stand  for  the  leading  terms  in  their  own  expansions 
for  0  large.  The  same  convention  (including  p)  in  the  remaining  equations 
(4),  (5),  (7)  completes  the  system.  The  relevant  solution  of  equation  (46) 
behind  the  reaction  zone  is 

(48)  T  =  T*. 

Note  that  the  equations  are  essentially  multi-dimensional.  Non-planar  effect 
are  no  longer  removed  to  distances  0(e),  where  diffusions  is  negligible:  and 
the  problem  is  correspondingly  more  difficult.  Note  also  that  the  actual  dist-jr.c  >:: 

over  which  changes  take  place  are  very  small,  being  measured  on  the  scale  of 
preheat  zone  (typically  less  than  0.1  mm.). 

To  complete  the  formulation  we  need  jump  conditions  across  the  reaction 
zone.  While  some  of  these  could  be  deduced  from  Sec.  II. 4,  the  generalisation 
in  notation  here  makes  it  more  convenient  to  obtain  all  of  them  ab  initio. 
Introduce  the  coordinate  (10)  at  the  point  of  interest  and  set 

(49)  n  =  e-1£. 

The  temperature  within  the  flame  zone  is  now  written 
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(50)  t  =  t#  +  e"1t1(5,y,z,t)  +  o(e"2) 

so  that 

(51)  p  =  -e  "1T;1t1  +  0(e~2)]; 

likewise  we  set 

.(52)  y  =  e“1y1(?,y,z,t)  +  o(e“2), 

(53)  v  =  y#(y,z,t)  +  e“1Y1U,y,z,t)  +  0(0~2), 

and 

(51*)  p  =  P0(S,y,z,t)  +  0{6~1). 

As  the  arguments  of  the  various  functions  imply,  the  temperature,  density, 
mass  fraction  and  velocity  are  continuous  to  leading  order  across  the  flame 
sheet  but  the  pressure  jumps.  The  goal  of  the  reaction-zone  analysis  is 
to  find  explicit  expressions  for  this  jump  and  those  in  the  normal  derivative 
of  the  ether  variables. 

The  variables  y^  and  t^  satisfy 

(55)  =  -32t1/d?2  =  (Y2/2T^)y;iexp(t1/T2) 

so  that,  in  particular,  y  +  t^  is  a  linear  function  of  £  which  matching 
shows  to  be  in  fact  independent  of  £  (to  leading  order  3H/3n  is  zero). 
Hence  h  is  continuous  across  the  flame  sheet  and  we  may  write 

(56)  yi +  ti  = 

In  the  usual  way,  integration  of  the  equation  for  t]_,nov  yields  sn 
expression  for  the  temperature  gradient  immediately  ah:  of 


the  reaction 
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zone  (coming  from  3t^/3£  as  £  -*■  -°°) ,  which  may  "be  written 
(57)  t|“)  =  -Y1  exp(h#/2T#). 

[  ]  denotes  the  jump,  i.e.  the  difference  "between  conditions  on  the  hot 

and  cold  sides  of  the  flame  sheet. 

We  turn  now  to  the  continuity  and  momentum  equations,  and  immediately  note 
that  to  leading  order  the  former  confirms  the  continuity  of  across  the 

flame  sheet  while  the  latter  confirms  that  of  Yl"  To  next  order  the 
continuity  equation  yields 

(V  +  v  3t  .  3v 

iv-0 

-1* 

so  that 


(58) 


3w  V  +  v  „ 

[  _£]  = - £_  [— )• 

1  3nJ  T*  3ri 


Similarly,  the  normal  component  of  the  momentum  equation  yields 

2 


3P  , 


4P  3  Vnl 


so  that 
(59) 


0  '  -  H  *  3*  V 


r  l  W'r  ^vn. 

tpl  ■  rt1' 


Finally,  the  other  components  of  the  momentum  equation  yield 


fl-ii 

at2 


=  0 


so  that 


3~L 

[  “  1  *  o. 

1  3n 


(60) 
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The  jump  conditions  are  completed  by  noting  that 
p(—  +  v  -v  )h  =  (1  -xe-1)v2ii  +  xe"1v2T 
is  an  exact  equation  everywhere, so  that 
(6l> 

In  particular,  we  see  that  when  the  Lewis  number  is  exactly  one  (x  =  0)  the 
function  h  will  vanish  identically  in  the  absence  of  inhomogeneous  initial 
or  boundary  data. 

For  slowly  varying  flames  it  does  not  matter,  so  far  as  stretch  is 
concerned,  whether  the  flame  surface  is  taken  to  be  the  hydrodynamic  disccr.ti 
or  the  reaction  zone  since  Vj^  does  not  change  across  the  preheat  zone.  On 
the  other  hand,  flame  speed  differs  for  these  two  possibilities  since  v^ 
does  change  across  the  preheat  zone.  For  near-equidiffusional  flames  the 
stretch  in  different  also,  because  of  changes  in  .  Nevertheless,  it  is 
natural  to  select  the  flame  sheet  rather  than  the  hydrodynamic  discontinuity, 
which  is  consistent  with  slowly  varying  flames  for  stretch  but  not  for  speed. 
The  determinations  of  speed  would  be  consistent  if  calculated  relative  to  the 
burnt  gas,  but  that  would  flout  long-established  convention. 

6.  Reduction  to  Stefan  Problems. 

The  problem  described  in  Sec.  5,  for  the  first  time  in  full  generality, 
presents  formidable  mathematical  difficulties.  In  general  a  numerical 
treatment  would  he  needed,  involving  the  simultaneous  solution  of  six  first- 
second-order  partial  differential  equations  (some  nonlinear)  on  each  side  of 
(possibly)  moving  discontinuity  which  has  to  be  found.  Even  the  description 
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of  small  perturbations  of  a  plane  flame  (Sec.  III. 8)  demands  the  solution  of 
equations  with  complicated  variable  coefficients. 

These  difficulties  are  significantly  reduced  when  the  fluid  mechanics 
is  eliminated,  in  the  spirit  of  Fig.  2  (Buckmaster  1979a,  b).  Then  only 
equations  (46),  (47)  remain  with  p  and  y  prescribed,  the  first  as  a 
constant  and  the  second  as  a  solution  of  the  continuity  and  momentum,  equations. 
This  is  the  constant-density  approximation  of  Sec.  1.5,  leading  to  the  system 

(62,  63)  (™  +  v  •  V)T  =  V2T,  +  y  -V)h  =  V2h  +AV2T 

when  the  density  is  absorbed  by  the  distance  and  time  units.  The  accompanying 
jump  conditions  are 

(64)  [||]  =  -X[— ]  =  XYiexp(h*/2T2), 

h  and  T  being  continuous  across  the  flame  sheet.  An  attractive  feature  is 
that  the  field  equations  are  now  linear;  but  we  are  still  faced  with  a  spatially 
elliptic  free-boundary  problem,  even  for  the  simple  case  of  a  uniform  flow 

(65)  y  =  (U,0,0)  with  U  -  const. 

For  simplicity,  consider  steady  plane  flow.  If  U  is  large,  i.e.  the 
gas  speed  is  much  greater  than  the  flame  speed,  then  a  flame  sheet  would 
have  to  be  rlr.ost  horizontal  (with  slope  of  order  U  ^);  under  the  right 
circumstances ,  changes  in  the  x-direction  would  then  be  much  smaller  than 
those  in  the  y- direction.  An  example,  discussed  in  detail  in  Ch.  IX,  is  that 
of  a  flare  next  to  an  adiabatic  wall  parallel  to  the  x-axis.  Under  such 
circumstances  the  field  equations  become  parabolic;  for,  by  writing 


(66) 


x  =  ux 
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in  equations  (  ) ,  we  obtain 

(67,  68) 


2  2  2 
3T  3  T  3h  _  3  h  3  T 

“  3y2*  "  3x*  3y2 


in  the  limit  U  -*•  ».  Moreover,  the  normal  n  is  almost  vertical  so  that  the 
n-derivatives  in  the  jump  conditions  (6U)  may  be  replaced  by  y-derivatives . 

The  result  is  a  generalized  Stefan  problem,  readily  amenable  to  numerical 
analysis  (Buckmaster  1979a).  Clearly  the  argument  is  net  affected  in  substance 
if  U  is  a  function  of  y;  the  x-derivatives  in  equations  (67),  (68)  acquire 
coefficients  depending  on  y,  the  same  for  each.  They  will  form  the  basis  of 
the  discussion  of  shear  flows  presented  in  Ch.  X. 

A  similar  reduction  is  possible  for  straining  flows,  i.e.  flows  in  which 
the  speed  varies  significantly  along  streamlines  (Buckmaster  1979b).  When  the 
speed  is  large,  the  angle  between  the  flame  sheet  and  intersecting  streamline 
is  small.  Then  changes  along  the  streamlines  occur  over  a  characteristic 
length  in  the  flow  field,  whereas  changes  normal  to  them  occur  over  a  flame 
thickness  in  the  neighborhood  of  the  sheet.  If  the  ratio  of  these  lengths 
is  large,  a  parabolic  description  resrilts. 

To  see  this  in  detail,  consider  a  curvilinear  system  of  coordinates  (s,n) 
in  which  s  measures  distance  along  one  of  the  intersecting  streamlines  and  n 
distance  from  it.  The  (iivergenceless)  velocity  field  is  supposed  to  have  the  form 


(69) 


v  =  Uq(a,v)  with  a  =  s/U,  v  ~  n/U, 


so  that  U  characterises  both  the  magnitude  of  v  and  the  distance  over  which 
it  varies.  Letting  U  »  gives  a  distinguished  limit  in  which  wy  is  0(l). 
The  combustion  field  is  described  in  terms  of  the  variables  a  and  n 


so  that,  in  the  limit 
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Vs  "  UV  Vn  =  ~nq0  Vlth  *0  =  %(a’0) 


is  the  velocity  field  experienced  by  the  flame.  Equations  (1*6), (1*7)  therefore 
simplify  to 


(70) 

(71) 


<*0  is  -  to,T  •  7§. 

3n 


(<1o  i;  ~  n  = 


32h 

3n2 


+X 


These  are  the  basis  of  the  discussion  of  straining  flows  in  Ch.  X.  The 
previous  equations  (67),  (68)  are,  of  course,  obtained  on  setting  q^  =  1,. 


o  =  X,  n  =  y. 
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